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Abstract: We define a replica field theory describing finite dimensional site disordered 
spin systems by introducing the notion of grand canonical disorder, where the number of 
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predicts a spin glass phase. We discuss the replica symmetric phase at length, explicitly 
identifying the correlator that diverges at the spin glass transition. We also discuss the form 
of continuous replica symmetry breaking found just below the transition. Finally we show 
how an analysis of ferromagnetic ordering indicates a breakdown of the approximation. 
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1. Introduction 

We consider site disordered models in finite dimensions with the following Hamilto- 
nian: 

H -- 



\Y,J{Ti-r j )S i S j (1.1) 
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in which N Ising spins are fixed at random points, r^, and are subject to a deterministic 
potential J(r). For example, a positive J(r) decaying with distance describes a dilute 
ferromagnetic system. Antiferromagnetic systems can also be treated, and because there is 
no lattice in this picture there will be no antiferromagnetic ordering. We also have in mind 
oscillatory RKKY-type interactions which cause frustration, as in the antiferromagnetic 
case. The physics we therefore address is related to the two types of ordering that can occur; 
ferromagnetic and spin glass, and to the relation between them. We shall discover that the 
approximation we use to solve the model is not completely reliable for the ferromagnetic 
system, nevertheless we present the general development for ferromagnetic interactions 
indicating where sign differences arise for the antiferromagnetic case. The Hamiltonian 
( |1.1|) could also be interpreted as an infinite range Sherrington Kirkpatrick model in 
which the bond strengths are correlated and have been chosen from a distribution very 
different from those we usually consider. 

This Hamiltonian can describe experimental systems rather well. However, analytic 
studies of disordered spin systems are usually based on lattice models in which the bonds 
take random values, such as the Edward Anderson Hamiltonian ||. Analytic work based 
on (|1 . 1|) has been hampered by the lack of a suitable field theoretic model (however a lattice 



based formulation has been proposed || ) . By considering a situation in which the number 
of spins in the system is random but quenched we are able to write a replica field theory 
for these site-disordered systems ||]. This field theory seems to be simpler than many of 
those coming from bond-disordered and diluted lattice models and should be accessible to 
many standard analytical techniques. 

In this paper, to gain an overall picture of the model, we use the technique variously 
known as Gaussian variational, Hartree Fock, Random Phase approximation and other 
names. It is necessary to use a method more sophisticated than mean field theory because 
the field theory that describes the model is expressed in terms of the magnetic variables 
and to understand the spin glass physics one has to look at composite operators. Mean 
field theory alone misses this spin glass physics, whereas the Gaussian variational method 
indicates that a spin glass transition is ubiquitous. We find that using the Gaussian 
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variational approximation the model can be solved in considerable detail, analytically in the 
high temperature region, and with the help of numerics at low temperature. The general 
picture that emerges is in accord with physical expectations, however there are some points 
that throw the reliability of the approximation into doubt. This is not surprising since 
although the approximation is expected to be good in the case of Heisenberg spins with 
many components, here we apply it to one-component Ising spins. The infinite component 
spin case is discussed in a companion paper ||. As we shall show, the problems with the 
approximation appear most clearly in two issues. Firstly, we shall find that for a purely 
ferromagnetic interaction a spin glass transition is predicted to occur at slightly higher 
temperature than the ferromagnetic transition and we present a proof of the impossibility 
of such a situation. Secondly, the approximation yields results that are not as dependent 
on the range of interaction or dimension as one would expect. This should act as a 
caution given the rather widespread application of the method in the literature. For this 
reason and also because a thorough understanding of this simplest non-trivial approach 
is a prerequisite for further advances, we feel that this study is worthwhile despite the 
deficiencies of the approximation. 

One of the interesting aspects of this work is that we are able to investigate the 
spin glass transition in a finite dimensional model. In the formalism we present, replicas 
enter the theory at an early stage and the spin glass ordering is intimately connected 
with the symmetry in replica space. The role of replica symmetry breaking (RSB) is of 
great interest in spin glasses, where the correct description of three dimensional materials 
is still controversial. Although the RSB in the mean field theory for spin glasses is now 
well understood || and related to the proliferation of pure states of the system, in finite 
dimensions the picture is less clear. Alternative qualitative approaches based on droplets 
view the spin glass phase as a disguised ferromagnetic phase with only two underlying 
fundamental states. The role of RSB in systems which undergo a ferromagnetic phase 
transition is maybe even less clear. We might mention the random field Ising model ||, 
dilute ferromagnets |J and the whole issue of renormalisation flow in the presence of RSB 
||1C|| . Experimentally, the problem is reflected in difficulty with the so-called reentrant 
transitions Jl l| . Part of the theoretical difficulty in these cases lies in the lack of a simple 
model in which the ferromagnetic transition can be explicitly analysed in the presence of 
RSB. All these issues provide strong motivation to study the model Ql.l| ). 
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Some of this work has already been briefly reported in [ffl, here we attempt a more 
pedagogic approach, discussing the basic issues more throughly and giving detailed ex- 
amples. The paper is organised as follows, Section 2 explains "grand canonical disorder" 
constructs the field theory and discusses some of the measurement issues. In section 3 
we perform the Gaussian variational analysis of the theory to derive the variational equa- 
tions. These equations are discussed for a general potential in sections 4 and 5 which 
concern the replica symmetric (RS) and RSB cases respectively. In the remaining sections 
we illustrate our results with reference to some particularly simple interactions, namely 
Yukawa and RKKY-type. A discussion of problems with this approach along with some 
speculation, appears in the conclusion. Two appendices give details of the solution of the 
four index correlator and of the stability analysis. 



2. Grand Canonical Disorder 

Firstly consider a model where the number of spins N is fixed: iV spins Si are placed 
randomly at positions uniformly throughout a volume V. This type of disorder we refer 
to as canonical disorder, as the number of particles is the same for each realization of the 
disorder. The spins interact via a pairwise potential J depending only on the distance 
between the spins. The Hamiltonian is the one written down in the introduction, 

H = -^J(r i -r j )S i S j (2.1) 

ij 

We shall proceed in the derivation assuming that J(r) is positive, thus giving rise to ferro- 
magnetic interactions. Later, we shall also discuss purely antiferromagnetic interactions, 
and at stages in the development will point out the sign changes necessary for a negative 
potential. A Hubbard- Stratonovich transformation expresses the partition function as 

Z N = J2j Vct>[detJf3]-hxp j j ^J-^r-r')^') drdr' + ^0(^)5^ 

(2.2) 

Employing replicas, we average out the site-disorder by integrating over the positions 
using the flat measure: J v Yl dri. 

Z n = j V<p a ldetJ(3]-%exp(-^ J J ^ <j> a {r) J~\r - r')4> a (r') drdr' 

i \ (2.3) 

+ iVlog^ J Tr expQ2(j) a (r)S a ) dr 
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Where we have introduced the trace over single site spins S a as convenient notation rather 
than write explicit cosh(p a factors. A field theoretic analysis of the above theory is compli- 
cated by the presence of the log term in the action. We overcome this difficulty by making 
a physically desirable modification to the definition of the disorder. In general one might 
expect the system to have been taken from a much larger system with a mean concentra- 
tion of spins per unit volume, p. A suitably large subsystem of volume V will thus contain 
a number of spins iV which is random and Poisson distributed: P(N) = e~ pV ^ p ^} . This 
distribution must be used to weight the averaged free energy, 



-[3F~ = -/3^P(A0Ftv = ^P(N)logZ N 

N N 

= limVP(iV) [ ~?kz± ) = i im (E_z±) 
n^o^ \ n I n^o \ n J 



(2.4) 



so we are led to define the partition function, S n = ^ N P{N)Z N . By analogy with the 
statistical mechanics of pure systems, we shall call this type of disorder "grand canonical 
disorder" . 

The theory is defined by the partition function for grand canonical disorder: 



E n _ e -pV 



V<f> a [detJ0\-^exp(--^ j j ^ <f> a (r) J~\r - r')(j) a {r') drdr' 



+ „/lY«p(l>(r)S.)* 



(2.5) 



Expanding the trace one sees that the leading term mixing replicas corresponds to a 
random temperature or random mass, familiar from bond disordered approaches, and that 
depending on the choice of interaction one might expect similar renormalisation group 
results |T^,[|. The simplicity of this form is evident, but the physical content is not 
immediately obvious so we now turn to a discussion of this point. 

2.1. Physical Operators 

In order to relate this theory to measurable quantities we must understand the type 
of excitations occurring in the theory and identify them with physical operators. This is 
not a simple task since the bound states of (|2.5|) are certainly not apparent at first glance. 



We start by returning to the original formulation of the model and considering the 
operator most closely associated with the field <j) a appearing in the theory. This is the spin 
density operator, 

M a {r) = Y J ^-r l )S? (2.6) 

i 

which appears in the replicated but unaveraged action as a source for (j) a . 

Z n = J V<t> a [detJp}-%exp ^--^ J j ^0 o J"V« drdr' + J ^M a </> a dr^j (2.7) 

The equations of motion (or Ward identities) at this unaveraged level provide relations 
between expectation values of M a 's and a 's. The equations of motion constitute a pow- 
erful tool and continue to hold for the fully averaged theory. For example the simplest 
relation is (M a ) = (<j) a ) / f3J(0). To make a connection with measurements performed on 
a particular sample with given disorder, we must look at self averaging quantities. One 
such quantity is the magnetisation density M, and below we show how it is related to the 
averaged theory, the replicated theory and finally to the expectation values of fields 4> a . 

M I /*(MM = P [<S>] ^ = lim ^ E / *<".(•■» = jjgj- (2-8) 

The over-line denotes the average over the full grand canonical disorder. In the third 
expression we have inserted the definition of the magnetisation density operator to show 
the relation to the spins, and in this case the square brackets indicate the average over site 
disorder alone. In the last term we have assumed that (<f> a ) is independent of the replica 
index. 

Similar arguments can be applied to the correlation function of magnetisation density 
operators. Firstly, (M a {r)M (r')) , is related to the field correlator (0 a (r)0(,(r')), by, 

Where we have assumed spatial translation invariance. In an experiment, the neutron 
elastic scattering cross section is proportional to the following correlator, 

(M(k)M(-k)) = p 2 \(S(k)S(-k))} = lim - Y^(M a (k)M a (-k)) (2.10) 

a 

The connected and disconnected parts of this correlation function can also be determined 
separately. 
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The spin density operator we have considered so far is a simple generalisation within 
the grand canonical context of the familiar spin operator. The expectation value, M in 
( |2.8|) , is the order parameter for ferromagnetism, but M a (r) is certainly not the operator 
sensitive to spin glass ordering. To probe this aspect of the physics it is natural to consider 
another operator: 

q ab (r) = Y^S(r-r i )StSl (2.11) 

i 

This operator is familiar from spin glass physics, its expectation value is a order parameter 
for the transition and its correlator is related to the non-linear susceptibility. In fact, bond 
disordered models such as the Edwards Anderson model lead to field theories [|13] in which 



the basic field variable is q a b{f) itself. These models immediately lead to non-trivial mean 
field theories describing the spin glass order parameter. The theory ( |2.5| ), on the other 
hand, contains only the simple fields (p a so the spin glass physics is hidden in the composite 
operators ( p.ll|) which do not appear at mean field level and necessarily require a deeper 
understanding of the field theory. In general we cannot provide such a simple and general 
relation between expectation values of q a b{f) and of the fields <p a as we were able to for the 
spin densities. However, within the context of the variational approximation, we will see 
that using linear response we can derive expressions for (q a b(f)) and equations obeyed by 
the correlators (q a b(k)q c d(—k)) . The connection between (q a b(r)) and the usual physical 
order parameter of Edwards and Anderson is clear from the relation (q a b(^)) = p[(^) 2 ] av - 
Although physical observables do not directly measure the correlators, these functions 
form a central part of our understanding in finite dimensions of both bond-disordered spin 
glasses and the site disordered systems considered here. We shall concentrate our attention 
on the connected part of the spin correlator which is related to the q a b correlators as follows: 



P 2 



(S(k)S(-k)) 2 con \ = lim 1 -^ ((foWtoHO) - (q aa (k)qbb(-k))) . (2.12) 

dytb 



More detailed information about the theory could be deduced by studying higher moments 
of the spin probability distribution. To do this, operators involving more spins can be 
introduced in a similar way. 
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3. Variational Method 

As we have already argued, mean field theory does not describe the interesting physics 
of this model and one must consider some improvement. In this section we analyse the 
theory using the Gaussian variational method which will allow us to substantially develop 
the theory without making an explicit choice for the form of the interaction J(r). The 
Gaussian variational equations are truncations of the full Schwinger Dyson equations and 
become exact in the limit of many spin components (such an m-component theory is treated 
in a separate publication ||). We find it convenient to phrase the development in terms 
of a variational method which allows us to calculate thermodynamic quantities simply, 
but from the outset, we would like to acknowledge the difficulties of variational methods 
for field theory [|H| . In the context of disordered systems, this method has had success 



in calculating exponents for random manifolds [15 1, and earlier in the problem of random 



heteropolymers [jTB . On the other hand, it has been less successful in situations where the 



RSB occurs in a one step pattern. One should also bear in mind that important effects 
may occur at higher orders in 1/m and indeed we shall discover that the approximation 
does not obey certain requirements expected on general grounds for purely ferromagnetic 
interactions. 

From the variational point of view, one selects a trial Hamiltonian, H t , and the method 
may be motivated as follows: 



J (3.1) 

> e -((H-H t )) t e -F t e - P V 

So the variational free energy, n(3F var = F tria i + (H — H tria i) tria i + pV , provides a bound to 
the true free energy: F% < F var . When the replica limit is taken, this bound is no longer 
rigorous, nevertheless the expression for the variational free energy is still valid. The 
Gaussian variational method simply consists in making a Gaussian trial, that is, choosing 
a trial Hamiltonian that is quadratic. Mean field theory can also be viewed as variational, 
in which case it corresponds to a linear trial Hamiltonian, and in this sense the Gaussian 
variational method is a simple generalisation of mean field theory. 

We allow the possibility of ferromagnetic order and make the following Gaussian 
anzatz (in which we have assumed translational invariance), 

Htriai = \j J>*( r ) - ^) G ab( r ~ r')(<t>b(r') ~ drdr' . (3.2) 

ab 



The variational parameters, (f> a and G a b(r) are simply related to <f) expectation values: 
((f) a ) = (f> a and ((f> a (r)(f> b (r')) = G ab (r — r'). This gives rise to a variational free energy of 
the following form, 



n(3F var = -^TrlogG ab + ^ J J ^ (0 a J~\r - r')(f> a + G aa (r - r') J~\r - r')) drdr' 
- pTr [e W (j2 $" S " + G - b ^ r ^ S - Sb ) dr + \ ^ lQ g(^ 3 ~ X ) + 

^ \ a ab / 

F /■ ^ r V f d d k n ,= ln , 



a 



(3.3) 



We have kept all the constant terms and have defined, 

= Tre H '= Yl ^p[^^aSa + l^2G ab (0)S a S b ) . (3.4) 

S a =±l V a ab / 

The variational equations follow immediately by varying with respect to (f> a and G ab . Their 
general form is: 

a J- 1 (O) = pf3n a = pf3Tr S a e H ' 

r-i ^ f- 1 no (3 - 5) 

where we have introduced generalisations of (|3.4|) : 



ms a s b 



n a = TvS a e H '= S a explj2^ S - + lj2 G - b 

S a =±l V a ab 

Q ab = TtS a S b e H = ^2 SqSb exp ( 4> a S a + ^ G ab (0)S a S b j 

S tt =±l V a ab / 



(3.6) 



Further generalisations involving more indices will appear in future equations. 

The variational equations, ( |3.5| ), can be interpreted diagrammatically as the leading 
equations of Schwinger Dyson, truncated in at the level of the four point function. As 
usual in this type of calculation, we find that Q ab is momentum independent; in many 
applications this only leads to mass renormalisation, but due to the replica limit non- 



trivial effects can occur \Wt. There exist some similarities with the formalism of 101. The 



next few sections will be devoted to solving these equations assuming various forms for the 
replica structure. 

The above derivation is for a ferromagnetic potential, in general the interaction poten- 
tial is not positive definite and to deal with it correctly one must split up the k space range 
accordingly. However there is a simple case, that of a purely negative or antiferromagnetic 
interaction. There will be no magnetisation, </> = in this case, and at the level of the 
free energy, the effect is a sign change of G ab in the definition of O (|3.4|). This leads to a 



reversed sign in the second variational equation and to a similar redefinition of fl a b- 

Using standard thermodynamic relationships we can also determine the entropy and 
internal energy from the free energy (|3.3|) , 



- + 5/ (2 -log/?J) (3.7) 

At high temperature, in the RS phase with zero magnetisation, these formulae simplify 
considerably. We find, S = pV log 2 as expected for Ising spins, and U = — pJ(0)/2 because 
in this limit the only spin correlation is from the same spin term of ( |2.1|) . 

3.1. Physical Operators in the Gaussian Approximation 

In the context of this approximation we can calculate expectation values of products 
of the physical operators identified in section 2.1. In general for variational theories these 
quantities should be determined using linear response as this leads, in the variational sense, 
to a smaller error than a direct evaluation [17]. At the level of ( |2.7|) , we introduce 



suitable sources for the operators of interest, follow through the analysis to obtain the 
generalised variational free energy, and then take appropriate derivatives before setting 
the sources to zero. Sources for the operators M a (r) and q a b( r ), only modify the free 
energy ( |3.3| ) through fi, which becomes, 

O(r) = Tr e H '^ = Tr exp ( £)(0 o (r) + (3h a (r))S a + \ ^(G a6 (r, r) + 0j ab (r)S a S b ) ) . 

\ a ab / 



(3.8) 



where h a (r) and j a b( r ) are the respective sources for M a {r) and q a b( r )- Translational 
invariance is lost until the sources are set to zero at the end of the calculation. 

To illustrate the method, consider the spin density operator M a (r). The magnetisation 
is defined as: 

M >>S^r A(r) - (3 - 9) 

Using the variational equation of motion for (p ( |3.5| ) we recover the equation M a {r) = 
\ J _1 </> a that we used in (|2.8|) . In fact, this method merely reproduces the equations of 



motion we had earlier for expectation values of spin density operators, and the correlation 
function is given by (|2.9| ). 

For the spin glass operator, q a b{ r )i there is no alternative method and in this case 
applying linear response gives: 

1 5F var 
P ojab{r) 



Qabcd(r - r') = (q a b(r)q c d(r')) 



^r, (3-10) 



1 S 2 F var p Sn ab (r) 



P 2 5j ab (r)5j cd (r>) (3 5j cd (r>) 

Using the variational equation of motion, some appropriate inverses and setting magneti- 
sations to zero we finally obtain the equations: 

Qabcd(k) = pQabcd + - ^ S 

abgh 

(k)Q ghcd {k) 

gh 

f d d p ~ ~ ( 3 - n ) 

Zabgh(k) =J2 n abef / T^W G eg (p) Gf h (k ~ p) 

ef 

Where £l a bcd is a trace of the form ( |3.4|) containing four spins. This equation goes beyond 
the simple Gaussian approximation, as is clear from its diagrammatic interpretation. It 
may be helpful for intuition to note that a similar diagram appears in the BCS theory of 
superconductivity. One reason for the importance of being able to get explicit expressions 
for the spin glass correlators is their use in showing the breakdown of the approximation 
for ferromagnetic interactions in section (4.2). 

Equations of this type, involving objects with four replica indices are familiar in the 
theory of spin glasses. In this particular case, note that although symmetry in a <-> b and 
c <-> d is manifest, the symmetry between pairs (ab) (cd) is not. Also remember that, 
in contrast to the SK mean field situation, diagonal terms with two equal indices do not 
vanish. In appendix A we solve this equation in the RS case but in the case of continuous 



RSB an extension of the methods of [USA would be needed. 
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4. Replica Symmetric Case 

We now specialise to Replica Symmetric (RS) solutions of the equations. Our param- 
eterisation of G a b is: go + 9i on the diagonal and, g± elsewhere. The constant term <fi a is 
taken to be replica independent. In this RS case we may write an explicit expression for Q 
( |3.4| ) , which leads to the following form (with ferromagnetic signs) for the variational free 
energy, 

F var 1 f d d k f _ 5 X \ 1 f d d k ~ , 1 T 2 f-i, n , 

" 2 9o{0) " J d ^ lQg ( 2 cosh(0 + ^v^(O))) + i J ^-(log/3J-l) 

(4.1) 

from which the RS version of the variational equations may be read off. 

According to the first of equations (|3.10|) , we define the Edward Anderson order 
parameter q as f2 & for a ^ b. This order parameter takes values in [0, 1] and is given 
by: 

q = n a ^ b = -L / ^e-^tanh 2 (^ + t^/g~Joj) (4.2) 
The variational equations for the components of the correlator are simply solved to give, 

9o(ft = = 

l-(l-q) P PJ(k) 

~m -V^ W/^(*0 (43) 



(l - (1 - q)pPJ{k)) 



2 



The value of the element g\ at zero spatial distance, <7i(0), appears throughout the equa- 
tions and can be determined from the momentum space form above as, 

f d d k P 2 J 2 (k) 

^ W (i _ (i _ g ) p 0j(A.)) 

The equation for is simply re-expressed in terms of the magnetisation M = <p/pj(0) 
( |2.8|) . Together with the expression for g, we obtain a pair of equations defining the two 
order parameters and which bear a striking resemblance to the RS mean field equations 
for the Sherrington Kirkpatrick model ffl, 



(4.5) 



M = -j= [ d£e~^ tanh ((3J(0)M + Zy/g^O) 
\Jlix J ^ 

q= [ ^e-^tanh 2 (p J(0)M + £y/g~i(ti) 
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In general a numerical solution is necessary, but at high temperature and low density a 
unique q = 0, M = paramagnetic solution is expected. Depending on the interaction 
potential, which may or may not allow a ferromagnetic state, the solution leads to two 
types of critical line corresponding to the order parameters M of ferromagnetism and to q 
of spin glass order which divide the phase diagram in the temperature density plane. These 
solutions however, may be either unstable or be energetically unfavourable with respect to 
solutions in which replica symmetry is broken. To analyse these issues in detail requires a 
choice of potential and a numerical solution, which we give for certain examples in sections 
6 and 7. Meanwhile we discuss the general conclusions that can be drawn from the simple 
q = 0, M = solutions, the stability criteria and issues that arise if there is ferromagnetic 
order. 

Given a solution of ( f4.5| ) for q and M, the correlators for the spin density operators 
can simply be determined by (|2.9| ) using the variational form: (<p a (r)(j)b(r')) = G a b- The 
correlator for the q a bS is less immediate and in Appendix A we give the general RS solution 

of (ED- 

4-1. Region with q = 0, M = 0. 



The M equation of (O) can always be satisfied by setting M = 0, and if a solution of 
the q equation exists then there is a paramagnetic phase. A particularly simple situation 
occurs in the region of high temperature and low density where q = is a unique solution 
for a wide class of potentials and dimensions. 

In this case, g± vanishes and the <p a propagator is diagonal in replica indices. It 
follows from (|2.9| ) that the magnetisation correlator is also diagonal. By using the results 

of appendix A for q = 0, or by noticing that for diagonal G a b, the second equation ( |3.11| ) 

d d P 

(2tt)° 

P 



simplifies to give T, abcd = Q a bcd J -4^d9o(p)9o(k - p), the q ab correlator can also be found, 



(M a {k)M a {-k)) 
(q a b(k)q a b(-k)} 



1 - p(3J(k) 



p ( 4 - 6 ) 



1 ~pJ ]Syi9o(p)go{k -p) 



All other q ab correlators vanish except for (q a a(k)qbb(~ k)) = p. The spin correlators follow 
simply from these formulae according to (|2.10| ) and (|2.12| ). 



Both correlators may, depending on the type of interaction, diverge somewhere in the 
temperature density plane. As temperature is lowered, the divergence in [(S(k)S(—k)) 2 ~\ av 
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will always occur before (or at the same point) as the divergence in [(S(k)S(— fc))] . If we 
assume that the first divergence occurs at zero momentum, and can therefore correspond 
to a continuous spin glass phase transition, the line on which it occurs is given by the 
condition: 

P J ^M k )Uk) = | = 1 (4.7) 

Where we have used the relation g\ = pqg^ Ql.3| ), which shows that g\ is linear in q. 
The spin density correlator generally remains finite on this line. The above condition is 
identical to the condition for the first appearance of a non-zero q solution in the equations 
Q4.5Q , since by expanding the right hand side of the q equation (for M — 0) we have, 

q = g 1 -2g 2 1 + %g 3 1 ... (4.8) 

The coincidence of the two lines is not surprising as it only relates to two different ways 
of identifying the line of the spin glass transition. 

In Appendix B we present the usual stability analysis along the lines of the work by 
de Almeida and Thouless , and find the following condition that the replicon eigenvalue 
remain positive, 

f d d k n r , . 

pr J Wf 9o = i gi<1 (4 ' 9) 

Where r is given by 

If ^ 

r = (1 - 20 bc + Q abcd ) = — / ^e-^cosh- 4 (/3J(0)M + ^v^(0)) (4-10) 

yzn J 

The solution discussed above, in which q vanishes, has r = 1, and is only stable above 
the same spin glass transition line, now also identified as the AT line. Just below the 
AT line, the perturbative RS solution with q small is always unstable. Performing the 
same expansion as in ( |4.8[ ) for small g± with M = for the r equation above, we find 
r = 1 — 2gi + 1g\ +0(gf) . The condition ( |4.9p is therefore always violated and the AT line 
signals a non-trivial breaking of replica symmetry. This indicates that we have a situation 
similar to that found in the SK model: namely a high temperature q = phase separated 
from a spin glass phase with non trivial RSB. Usually the coincidence of lines noted here 
indicates that the RSB will take place continuously. The reader will doubtless note that the 
above analysis seems remarkably general both in terms of the type of interaction involved 
and on the dimensionality of the space. One sees that when presented with the statistics 
of the matrix J(n — rj) it is very difficult to divine from what interaction type and what 
dimension of space it arises from. The consideration of a genuine cut off interaction where 
one would expect the geometry and type of interaction to play a much stronger role is 
relegated to a future study. 
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4-2. Magnetisation 

For a ferromagnetic interaction the mean field equations, M = ptanh(/3J(0)M) sug- 
gest that we should expect a ferromagnetic transition at p(3J(0) = 1. The Gaussian 
variational approximation sees this transition as the divergence in the M a correlator in 
QL6|) . However it is clear from the preceding discussion that the spin glass transition in- 
tervenes and that we should not even be attempting to solve the RS equations ( |4.5|) to 
locate a ferromagnetic transition because the system is already in an RSB phase. The AT 
line shrouds the critical region and the RS equations are not relevant to the transition. 

Such a phenomenon may be surprising and the order in which the transitions occur 
deserves further discussion. For a generic interaction and generic number of spin compo- 
nents it is reasonable for the spin glass transition to occur at higher temperature than the 
ferromagnetic one. For example this is known to happen in the random field Ising model, 
where replica symmetry breaking has been shown to take place before the ferromagnetic 
transition ||. The problem occurs for Ising spins and a purely ferromagnetic interaction 
(positive in space) where we might well expect the transitions to be concurrent. 

Let us try and clarify this point [pOj by noting that the divergence in [(S(k)S(—k)) 2 ~\ ^ 



occurs while [(S(k)S(—k))\ is still finite. The divergence implies that at sufficiently 
large spatial separation the spin glass correlator is larger than the ferromagnetic one: 
[(S(r)S(0)) 2 ] av > [{S(r)S(0))\ av We might expect that these correlators are distributed 
according to some probability distribution P(C) and therefore write, 



(S(r)S(O))] = / P(C)CdC 

Jav J 

(S(r)S(O)} 2 } = [ P(C)C 2 dC 

J av J 



(4.11) 



But for Ising spins we expect that the correlator C = (S(r)S(O)) is bounded, < C < 1, 
provided all the interactions are ferromagnetic. Using such bounds in the expressions 
above, it is clear that the inequality, [(S(r)S(O)) 2 ] av < [{S(r)S(0))~\ holds for all r. 
This inequality is violated in our solutions. It should be possible to trace this breakdown 
of the Gaussian variational approximation to neglect of some higher order (in number 
of spin components) term in the Schwinger Dyson equations. Indeed, Sherrington [pT 



has pointed out the dangers inherent in methods similar to the one we use in the case 
of a Landau Ginzburg approach to dilute ferromagnets, and has identified the relevant 
diagrams. 
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Even though it is inappropriate to look for a continuous ferromagnetic transition in 
the RS equations, it is interesting to see what happens. The usual argument, based on the 
graphical solution of (|4.5|) , provides a simple condition on the gradient of the magnetisation 
equation (at M = 0) for the existence of small M solutions: (1 — q)p[3J(0) > 1 (the 
(1 — q) term causes a suppression of the transition temperature with respect to the mean 
field prediction). On the other hand, for the Gaussian integrals to be well defined the 
propagators must be positive leading to, (1 — q)p(3J(k) < 1, V7c. In this sense the usual 
condition always fails, even in finite volume. This argument does not preclude RS solutions 
with magnetisation, it only indicates that the ferromagnetic transition itself is hidden, or 
possibly first order. Indeed for suitable ferromagnetic interactions numerical analysis of the 
equations allows us to find ferromagnetic solutions below the transition with no difficulty. 
At zero temperature general arguments would require such a solution for interactions that 
do not vanish beyond some range. Since this issue requires a choice of potential we will 
discuss the issue further when we look at examples. 

5. Replica Symmetric Breaking 

Since we have found a scenario similar to that of the Sherrington Kirkpatrick model, in 
which the RS q ^ solution is unstable as soon as it appears, we shall look for continuous 
replica symmetry broken solutions. We consider Parisi matrices in which the off-diagonal 
part of the matrix G ao (k) is parameterised by a continuous function g(u,k) where u G 
[0, 1], and the diagonal part is denoted by #£>(&;). The algebra of such matrices has been 
developed, and expressions for products, inverses etc are given in the appendix of reference 
||15|| . We shall follow the notation used in that paper for certain integral transforms that 
occur frequently. Although the similarity to the random manifold problem seems clear 
the sign of the potential is opposite and is related to a random manifold problem with 
imaginary noise. 

For G a b a Parisi matrix, O ( |3.4| ), is very similar to the free energy in the SK model. 
It is well known that this cannot be obtained in a closed form and a standard strategy 
is to work close to the transition line by expanding O up to a term of 0(g 4 ) which in 
the SK model is the first term leading to a breaking of replica symmetry ||22|| . We shall 
assume that there is no magnetisation in this calculation because the denominators of the 
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magnetisation correlator showed no tendency to diverge at the AT line. The expansion is, 
Q 1 ) 9 n(0) - \ fig 1 + V - V - g [ U 9 2 ) du (5.1) 



n T v ' 4 Jo V 6" 3' 
where g is the spatial function, g(u, 0), evaluated at zero distance. The remaining terms 
in the action are easily computed within the algebra of Parisi matrices. The variational 
equations are, 



[g D (k)]- 1 = +pa D (k)=p[^-l 



(3p 



(50 

[g(u, k)}' 1 = -pa(u) = 2p— 



(5.2) 



Where the notation indicates the component of the inverse matrix which is also of 

Parisi form. In deriving these equations of motion it is important to take care of the signs; 
for example a minus sign appears in 5/5g(u, k) TrlogG = — [g(u, fe)] -1 , and the signs for 
a have been chosen in order that the following equations for the inverses appear in the 
simplest form. For the antiferromagnetic sign interaction the definition (|5.1|) is changed 
and a sign appears in the second of the variational equations. 
Defining the denominator combinations: 

D D (k) = a D + (a), 

(5.3) 

D{u,k) = a D + (a) + [a](u), 

where following the notation of [|15| the angle and square brackets denote (a) = j^dua(u) 
and [a] (it) = ua(u) — f Q dva(v). The equations can be inverted to find, 

~ m _J__ f, f X du [a}(u) , ff (0) 



pD D (k) V Jo u 2 D(u,k) D D (k), 



pD D (k) \uD(u,k) J v 2 D(v,k) D D (k) 

Taking the derivative with respect to u of the second of these equations leads to the simple 
relation, 

8-=!' ±£1 (5.5) 
p D z up du D 

where we have made use of D' = [a]' = ua' . Proceeding by differentiating the second 
equation of (|5.2|) with respect to u and using the expression for g' above, one obtains 
a' = 0, or, 

l + g 2 -ug- [ g] = p (fAJLl.) . (5.6) 
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Taking another derivative in some region where equation (5J3) holds we find: 

/ r d d k 1 \ 



g = a(u)u 



u 
2 



l + 2p z 



r d a k l 



\ 



(/ 



d d k 



(27r) d D 2 J / 



(5.7) 



0.2 0.4 0.6 0.8 



Fig. 1: Typical RSB pattern just below the AT line. 



Substituting this result into the relation ( p.5|) we find that a is determined by a a 
non-linear, but first order, differential equation in u, 



^-(ua(u)) = — - [ 
au V / up J 



A d 1 (5.8) 



(27r) d du D 



We cannot go further without being more explicit about the form of the interaction and 
dimension. In general the function at(u) has a complicated dependence on a(u) and the 
differential equation will be difficult to solve, but based on the examples we have studied, 
we expect that the RSB pattern as shown in fig. 1 to be generic. The initial part of g(u, 0) is 
determined by a power series analysis of the nonlinear differential equation about the origin, 
and in all examples we have studied the leading term is linear. Above some breakpoint, uq 
close to the origin, g(u, 0) becomes constant. For consistency with the original expansion 
of O we must work close to the AT line and must require that the g(u, 0) is perturbatively 
small in deviations from that line. Other patterns could be envisaged, though we have not 
yet found any examples. 



17 



For certain choices of interaction or dimension the differential equation may simplify 
radically. For example in large enough dimensions (the critical dimension depends on the 
large momentum behaviour of the interaction; if J ~ k~ 2 , then d cr u — 4) the second term 
in a will disappear in the limit in which the short distance cutoff is removed. The resulting 
equation is simple and we find a scenario very similar to that found in the SK model. Due 
to a remarkable cancelation a is also a constant for Yukawa interactions in 3 dimensions. 

This analysis near the AT line is insufficient for some purposes such as to investigate 
continuous RSB solutions with magnetisation, and one would like to solve the equations 
in greater generality. More terms could be kept in the expansion (|5.1| ) |fL3f ||23|| , or it may 



be possible to use the differential equation obeyed by O along the lines of the SK case [ 5|] . 
Further difficulties would occur if it became necessary to consider the order parameters 



with more replica indices that can be constructed as in section 2.1 p5 ]. We do not attempt 
such generalisations here. 

We have also studied the equations for one step of replica symmetry breaking. We do 
not expect such solutions to be relevant near the AT line, and at lower temperatures their 
importance should be judged relative to some continuous solution. Another use of such 
solutions is in calculating the entropy of the meta-stable states which is relevant to the 



dynamics. We have made a calculation following Monasson |2(J that demonstrates that 
as usual in SK-like situations, this entropy remains zero down to the AT line and that we 
therefore do not expect any distinct dynamical transition. 



m=± — 2 , (6.i) 



6. Yukawa Potential. 

A simple potential for illustrative purposes is, 

1 1 

where the signs refer to ferromagnetic and antiferromagnetic interactions respectively. This 
momentum space form describes a Yukawa type potential that is screened on a length scale 
1//U. We choose to measure distances in terms of this scale (for example the dimensionless 
density becomes p/fJ- d ), and therefore set fi = 1 in what follows. 

The ferromagnetic sign leads to a rich phase diagram and the complications discussed 
in section 4.2. For this reason it is useful to consider the antiferromagnetic sign since 
antiferromagnetic order cannot exist in the absence of a lattice, and the phase diagram 
will be simpler, allowing us to concentrate on the properties of the spin glass transition. 

We discuss different dimensions separately, starting in three dimensions. 
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6.1. Three dimensions 

In three dimensions the real space form of (|6.1| ) is the well known Yukawa form, 

J(r) = ±f-^- (6.2) 



In this case the integral ( |4.4| ) can be done to obtain gi 



(/l ''' 7 ' (27r) 3 (fc 2 + 1 T (1 — q)pf3) 2 87Ta/1 =F (1 — <?)p/5 1 ' J 



which allows us to numerically solve the equations (|4.5| ). At high temperature the unique 
solution is q = 0, M = 0, whereas at lower temperatures for the ferromagnet, the positivity 
of go restricts the possible range of q to, q > 1 — l/p/3. 

The spin glass transition line is defined by the condition 8n^l =p (3p = p(3 2 , which can 
be solved to give, 

p = 2(4itT) 2 (tT + VT 2 + 1/(4tt) 2 ) . (6.4) 

This is the only phase transition in the antiferromagnetic case whereas in the ferromagnet 
it lies very close to the line, p/3 = 1, where ferromagnetic ordering would occur in mean 
field approximation. 

Above this line, in the high temperature phase with q = 0, M = 0, only the diagonal 
magnetisation correlator remains: 

(M a (r)M a (0)) = pS(r) + 9 P ^ (6.5) 

We have given the real space form because it has a simple interpretation. The first term 
is from the same spin, and the leading piece of the second term arises from a two spin 
term. Only considering two spins, the correlator can be written in terms of Boltzmann 
factors which when expanded for small interaction strength give p 2 f3J{r). Using ( |6.2[ ) one 
recovers the leading piece of the full correlator. 

Still in the high temperature region, we can also calculate the only non-trivial q a b 
correlator for the P combination of indices (see Appendix A): 

Qabab(k) = (q a b{k)q a b(~k)) 



1 ~pf (£y^ o 2 (*0 ! _ p§1 i sin -i ( * ^ 

"' 4tt k y v /4(l T /3p)+2fc 2 J 



(6.6) 
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The connected spin correlation function [(S(k)S( — k)) 2 on ] av is obtained by subtracting the 
term (q aa {k)qbb(~k)) = p ( |2.12| ). In the small momentum limit this expression simplifies 
and becomes proportional to l/(k 2 + m 2 ) with mass given by, 

In accordance with the general analysis, criticality occurs for vanishing m 2 which is another 
way of finding the AT line. Despite the fact that we are employing a Hartree Fock method, 
which in pure models can give non-classical exponents, here m 2 is an analytic expression 
and the exponent v takes its mean-field value of |. 
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Fig. 2: The regions in the density temperature plane where various RS solutions 
exist. For ferromagnetic Yukawa interactions in 3 dimensions. A region between 
the spin glass and ferromagnetic transitions is unphysical. 

At lower temperatures, other replica symmetric solutions exist. The simplest case is 
the antiferromagnet for which there are two solutions below the AT line. One solution, 
q = 0, is unstable both with respect to longitudinal and replicon fluctuations, and the other 
has q > and is numerically found to be unstable with respect to replicon fluctuations 
as expected from the general analysis in section 4.2. The ferromagnetic case is more 
complicated, besides M = solutions (the q = solution only exists for (3p < 1) that 
are unstable just as for the antiferromagnet, there are solutions with M ^ 0. These 
magnetised solutions first appear on a second line at slightly lower temperature than the 
AT line, as shown in fig. 2. On this line there is one new stable solution with non- vanishing 
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magnetisation and q, these values being largest at the low temperature end of the curve. 
As the temperature is reduced the solution splits into two; the branch with larger values 
of M and q remains stable and the magnetisation continues to grow in the usual way, the 
other branch decreases in magnetisation and becomes unstable with respect to the replicon 
not far below the line shown on the figure. We suspect that the decreasing solution will 
not play any role. 

The significance of the upper line in fig. 2 and the new solutions that it heralds de- 
pends on whether continuous solutions also exist in this region. Free energies must then 
be compared to identify a first order transition. As was discussed in section 4.2, for purely 
ferromagnetic interactions we expect concurrent spin glass and ferromagnetic transition 
lines, and therefore we are witnessing a shortcoming of the approximation. It seems likely 
that at sufficiently low temperature the single stable new RS ferromagnetic solution will 
have the lowest free energy and will represent the system. If this is the case, at zero tem- 
perature the transition must occur for finite density, suggesting a connection between the 
ferromagnetic and a percolation transition. This would be a reasonable conclusion for a 
ferromagnetic potential that strictly vanishes beyond some radius (it would be interesting 
to consider such a potential to clarify this connection). On the other hand, for an inter- 
action of the type ( |6.2| ) which never vanishes, we should always expect a ferromagnetic 
RS ground state at zero temperature. Investigating this region in detail, we find that 
the upper line shown in fig. 2 does not quite extend down to zero temperature because 
the main magnetic solution becomes unstable with respect to replica fluctuations. This 
observation can be confirmed by making a low temperature expansion of the RS equations 
( |f .5|) . Apparently this is another failure of the approximation, and is maybe not surprising 
since it is unlikely that the (3 — > oo limit commutes with the large number of components 
limit used to justify the approximation. 

We now turn to the form of continuous RSB just below the spin glass transition. In 
view of the arguments showing that the critical region for the pure ferromagnet is not 
treated well by the approximation, we only discuss the antiferromagnet. 

It is useful to identify mass parameters by writing the denominators ( |5.^| ) as, 

D D (k) = j-(k 2 + m 2 D ) 

ml={l + (3p(l + (a))) 

D(u,k) = J- (k 2 + m 2 (u)) 
PP 

m 2 (u) — m 2 D + (3p[a](u) 
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Performing the momentum integrals over inverse powers of the denominator we find a 
remarkable cancelation so that the parameter a ( |5.7| ) is independent of u, 

1 (32n 2 
2\W 



« = o ( ~ZW ~ 1 ) ( 6 - 9 ) 



We consider solutions consisting of a continuous piece up to a breakpoint uq, followed 
by a constant as shown in fig. 1. The differential equation is trivially solved to give 
m(u) = rriD + 2irau 2 /f3. Continuing to solve the RSB equations we find that in the 
continuous region, a contains linear and cubic terms in u and that the breakpoint is 
determined by a quartic equation. Since uq is small we only give the leading terms, 

^(32tt 2 ± p(3 3 )u = (p/3 2 - 8n^l + rf?) . (6.10) 

The right hand side vanishes on the AT line, recovering the RS solution with uq = 0. 
Expanding to first order about this line we find an expression for the breakpoint in per- 
turbation theory, 

_ 5(3 0^)0^-160 
U °~ 4 J (p/33-32^ ^ 

This indicates that we have a consistent solution since g remains perturbatively small so 
the expansion near the AT line is justified. The magnetisation correlator can be calculated 
in this region along the lines of (|2.10| ), in which the sum now becomes an integral ( |5.4j) . 
To leading order we find 

a, (t)= ) (6.i2) 

k 1 + m z D \ p{k z + m A D ) J 

Where wid = VI + pP + O(uo). 

These detailed predictions of a spin glass phase for a site disordered antiferromagnet 
are interesting because few analytic results exist for the problem. Preliminary numerical 
simulations of the dynamics [^7| reach the same conclusion as those described by McLe- 
naghan and Sherrington |28| (who also refer to experimental realisations); that is to say, 



no such phase is observed. More work is needed in this area since for the analytic work, 
either the Gaussian variational approximation or the simple choice of a two replica order 
parameter may be in doubt; whereas the numerical simulations suffer from small sizes and 
consequent dependence on the boundary conditions. 
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6.2. Yukawa in d > 4. 

In four or more dimensions many of the integrals diverge indicating dependence on 
details of an unknown short-distance theory. For example, in the RS phase the integral 
defining <7i(0) diverges, the leading term (in the momentum cutoff A) gives: 

A-"- 1 



"=^(5=4) UJ <6 ' 13) 
Where T d = tt d /(2n) d = 2/ ((4n) d / 2 T(d/2)). The AT line is thus given by, 

So as the cutoff is removed (A — > oo), the AT line moves towards the p = axis, and the 
RSB phase covers the whole phase space. One can proceed to expand about this line to 
find a simple RSB with a = 1/2 and uq — 48(3/(3. 

Alternatively, we can recover a finite temperature transition by taking the infinite 
dimension limit while scaling the cutoff as A//i ~ \fd. This limit sets gi oc (3 2 qJ 2 and 
yields the SK mean field equations. 

6.3. Yukawa in d < 3. 

The two dimensional example contains most of the features seen in low dimensional 
Yukawa models, and we concentrate on it here. Again the integral can be done to obtain 
a simple form for gi, 

R 2 [ d 2 k 1 pH\ , v 

91 PPq J (2n) 2 (k 2 + 1 =p (1 - q)p(3) 2 4n (1 =F (1 — q) P 0) ^ ° J 

The AT condition is 4tt{1 =f p(3) = p(3 2 , so the line is given by: 

4nT 2 

» = T±^T (6 ' 16 » 

Note that in the antiferromagnetic case, the line always lies below T = 1/4-k. Numerical 
investigations of the RS equations lead to a picture qualitatively the same as the one 
described in the three dimensional case. Quantitatively there are differences, notably, in 
two dimensions the unphysical region between the AT line and the line on which stable 
ferromagnetic solutions first appear, is wider than in fig. 2. 
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^/-l / = ° ( 6 - 17 ) 



In searching for a RSB solution for the antiferromagnetic case it is useful to work in 
terms of the mass parameter m(u), defined by the combinations ( |6.8|) rather than directly 
in terms of a. In two dimensions we find a = !(^p-^ 2 — !)• The differential equation 
( |5.8|) can be written in terms of the variable x = u 2 , and it is possible to work in terms of 
/ = m 2 , to obtain, 

32tt 2 p\ df /16tt 2 
~W Xj ~n) dx~ + {~pJP 
Provided the breaking pattern is similar to that shown in fig. 1 we need only analyse the 
small u behaviour of this equation. We find results very similar to the three dimensional 
case: m? ~ m 2 D {l + nu 2 (16'K 2 m 2 ~ ) — /9/3 3 )/p/3 4 ), and g approximately linear in u. This 
is a consistent procedure since near the AT line we obtain a sensible expression for the 
breakpoint, 

"°= 4 ^!^!f- < 6 - 18 > 

Closer inspection of the differential equation shows that it is ill-behaved because so- 
lutions tend to blow up at finite values of the parameter u. Also note that two distinct, 
constant solutions exist, which may suggest a solution with one step of RSB. 

One dimension is interesting because exact results are known about the bond disor- 
dered version of the model. The AT line is given by 4(1 =f p/3) 3 / 2 = p/3 2 , which is a cubic 
equation in p. A series solution of the differential equation shows behaviour like that of the 
two dimensional case. The zero dimensional case is formally similar to an RKKY potential 
in any dimension, and is better discussed in that context. 

7. RKKY-like Potential 

The RKKY potential itself, which in three dimensions is J(r) = cos /xr/r 3 , has a 
complicated Fourier transform. So as is usual in analytic work we consider instead: 



J(fc) = Po/^, fork<p; (71) 
L 0, for k > /x. 



J(r) oscillates in sign to have the correct qualitative form, but note that the negative parts 
are not very strong. As was the case with the Yukawa potential, \x and Jo can be scaled 
away by measuring distances and temperatures appropriately, so we will set \x = Jq = 1. 
The integral for gi can immediately be performed, 

= 2 f d d k 1 _ T d p(3 2 q 

91 PP q j lkl<fl (2ny {1 - {1 - q )pp) 2 (l-(l- q ) p p) 2 

Where T d = tt d /(2n) d = 2/((4n) d / 2 T(d/2)), and in three dimensions, T 3 = 1/6tt 2 . 
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Fig. 3: Regions of existence of RS solutions for RKKY-like interactions in three 
dimensions. 



The AT condition is given by, 

T ^ 2 , = 1 (7.3) 

So at this transition the magnetisation propagators remain finite. This can be solved to 
give the AT line, 

P=(T+Y-^T+^f) ( 7 - 4 ) 

Since the form of this potential does not depend strongly on the dimension we shall 
discuss only the three dimensional case. Above the AT line, in the q = RS phase, the 
correlators are 

<*.(*)*.(-*)> = {j^ tltl; < 7 ' 5 ' 

and 

Q iab )(a b )(k) = { i-^(t^) (*+*K*-2) a ' I '' l ' < (7.6) 
L p, for k > 2. 

Where the non-trivial form of the second propagator arises from the slightly complicated 
integration region and this allows us to obtain a non-mean field exponent, rj = 1 . 

Numerical study of the RS equations exposes a scenario very similar to the ferromag- 
netic Yukawa case. The solutions with zero magnetisation are always unstable below the 
AT line, but new magnetised solutions appear on another line. We recognise the second 
line from the ferromagnetic experience, at very high densities the positive short range part 
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of the potential dominates and the system magnetises, returning to a RS form. The phase 
diagram is shown in fig. 3. 

For RSB we find a = |(1 + r >2 2 p/33 m 6 ). The differential equation can be written in 
terms of / = m 2 where m is defined in a similar way to 



T ^ xf - md y^ + ^ f ^ f=0 ( , 7) 

The expression for the breakpoint is sensible and given by, 

"° = iT ^ (i + jfd- P 0) ■ (T - 8) 

From this information we can calculate the structure of the two index correlators g(k : u), 
leading to non trivial momentum dependence in gr,{k), however the analysis only holds 
close to the spin glass transition and cannot address the ferromagnetic transition itself. 
Indeed, at the ferromagnetic line, g will be large and the expansion of O is manifestly 
inappropriate. In addition the four index correlation functions Q a bcd contain much of the 
physics of the spin glass phase: for example the ^-exponent may be extracted from the 
long distance behaviour of such objects. Equation (|3.11|) is however an equation carrying 
four replica indices and the solution in the case of continuous replica symmetry breaking 



is technically rather formidable requiring extensions of the methods described in 18 



8. Conclusions 

The Gaussian variational method, which is the simplest non-trivial way of analysing 
this site disordered field theory, provides an interesting picture of the model. For all forms 
of interaction between the spins we find a spin glass transition which separates a high 
temperature paramagnetic phase from a low temperature phase with non-vanishing Ed- 
wards Anderson order parameter. The transition is signaled by replica symmetry breaking, 
which takes place continuously just below the transition. We have carefully displayed the 
connection between physical observables and the replicated quantities we calculate and 
see explicitly how the spin glass correlator, (q a b(k)q c d(—k)) , diverges as the transition is 
approached from above. 

The properties of the low temperature region depend on the precise form of the in- 
teraction, but not sensitively on dimension in fewer than four dimensions. Our study of 
antiferromagnetic Yukawa interactions indicates that the low temperature phase breaks 
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replica symmetry continuously at all temperatures. The detailed form of the breaking 
has, however, only been determined in the vicinity of the transition. This prediction is 
interesting in view of the numerical simulations which contradict it and further work is 
needed in this case. In the case of ferromagnetic interactions, or the RKKY-like interac- 
tions which contain a substantial ferromagnetic component, we have found stable replica 
symmetric solutions with non-vanishing magnetisation. It seems likely that at sufficiently 
low temperature one such solution will describe the system. We are unable to be certain on 
this point until we know about the existence and properties of solutions with continuous 
replica symmetry breaking in the region away from the spin glass transition line. This 
problem also means that we are not presently able to directly analyse the ferromagnetic 
transition, though, because the magnetisation correlators remain finite at the spin glass 
transition, it is clear that the transitions are not coincident indicating a breakdown of the 
approximation in the case of a pure ferromagnet. Some further difficulties with the zero 
temperature limit have been mentioned in the text. 

Besides the inaccessibility of the ferromagnetic transition, an important aspect of the 
spin glass physics is missing. Namely, the behaviour of the spin glass correlator below the 
transition where it should remain massless yet have a non-trivial exponent 9 [|7j . We have 
derived equations obeyed by the correlator, but because it is a four index quantity, the 
technical difficulties of solving the equation in the continuously broken spin glass phase 
are presently beyond us. Both the problems mentioned can be addressed in the case of 
m-component Heisenberg spins in the limit m — > oo ||. 

The two issues in which the approximation is clearly unreliable are in its treatment of 
the purely ferromagnetic interaction and it relative insensitivity to dimension or interaction 
type. It therefore seems appropriate to add some further words on when we expect the 
approximation to be trustworthy. The Gaussian variational approximation can be seen as 
making a Gaussian ansatz on the Schwinger-Dyson equations of the theory for the fields 
4> a {x) ( m the case where there is RSB the anzatz is actually more refined in terms of nested 
Gaussians). However this ansatz can only be justified in certain circumstances. Some form 
of the central limit theorem must be brought into play, and the possible ways in which it 
would be justified are the following: 

• There are many effective neighbours for each spin giving a large sum contributing to 

the local field. This will be true in a number of cases such as large spatial dimension 

(where one is of course closer to mean field) and long range interactions for which even 
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in dilute systems there are many effective neighbours. Therefore one would expect 
that for certain short range applications the theory may be expected to fail. 
• When there are long correlations on the local field distribution one expects the central 
limit theorem and thus the approximation not to be applicable. Indeed for purely 
ferromagnetic interactions the correlation between local fields will become stronger as 
the temperature is reduced, we see the mass in the two-point correlator does become 
rather small in three dimensions even if it does not become zero before the predicted 
spin glass transition. In the non-ferromagnetic cases the correlations in the local fields 
do not become so strong and hence one may have more confidence in the method. 
It is natural to enquire what we could say about the system beyond the level of the 
Gaussian variational analysis. Simulations of the basic Hamiltonian we started with are at 
an early stage P7| , and the non-trivial dynamics usually associated with disordered systems 
has not yet been observed. The renormalisation group could yield important information 
about the model. A naive application addresses the ferromagnetic transition because these 
are the obvious variables in the problem, but a perturbative approach leads to exactly the 
problem investigated by Dotsenko, Harris, Sherrington and Stinchcombe ||. These authors 
found that the renormalisation flow in a dilute ferromagnet was disrupted into a replica 
symmetry broken form. The precise interpretation of this observation is as yet unclear, 
but we hope that the analysis we present in this paper based on the Gaussian variational 
method will shed light on the problem. One might hope that a more sophisticated approach 
to the renormalisation, in which the appropriate variables for the spin glass transition were 
isolated, would give new and interesting information. These issues we postpone to another 
publication |L2 . 



We have observed continuous RSB in finite dimensional models, and providing an 
interpretation by identifying the pure states in physical terms would be an interesting 
task. For this purpose, as mentioned earlier it will be interesting to consider a potential 
that strictly vanishes beyond some radius, as the connection with percolation is clearer 
and one would expect to see a greater sensitivity on spatial dimension. 
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Appendix A. The RS correlator for q a f, 

In this appendix we resolve the equations ( |3.11| ): 

Qabcd{k) = P^abcd + - ^ S 

abgh {k)Q g hcd{k) 

gh 



r d<i p _ _ (A- 1 ) 
^abgh(k) = 22tt abef / —— 1 G eg (p)G fh (k-p) 

ef J [Z7t) 

for the correlator Q a bcd(k) = (q a b(k)q c d(—k)) in the general paramagnetic RS case with 
M = 0. but arbitrary q. 

Notice that although symmetry in a <-> 6 and c <-> d is manifest, the symmetry between 
pairs (ab) <-> (cd) is not. This could be rectified by an appropriate matrix multiplication, 
but it is simpler to leave. 

These equations are a set of linear equations for the nine possible index combinations 
we have to consider since the diagonal terms are inevitably mixed into the equation (the 
missing symmetry would reduce this to the seven combinations, A, B : C, D, P, Q, R in the 
notation of jn|). The four spin trace Q a bcd is completely symmetric and takes values either 
1 or q except for the case with all indices different where it is equal to — l + 2q+r. T> a b g h(k) 
is defined by quadratics in g and g\. It is definitely not pair symmetric, but does satisfy 
S( aa )( aa ) = S(aa)(6b) and E( ao )( a 6) = £( aa )(6c)- Other simple combinations are based on 
the longitudinal, (P — 4Q + 3R), and replicon, (P — 2Q + R), mixtures of P: S( a 6)( a 6) 7 Q- 
S( a 6)(ac) and R: S( a 6)(& c )- 

The 9x9 matrix equation splits into a 4 x 4 and a 5 x 5 block. The 4x4 part is 
simple giving, Q( aa )(aa) = Q(aa)(bb) = P- The 5x5 part contains the P, Q and R pieces 
and the longitudinal and replicon combinations of the Q's are particularly simple, 



Qlon = Pi 2 + 2g + 3r) , D lon (k) = l- P (-2 + 2q + 3r) [ ^~g (k - p)(g - 2~ 9l ){p) 

Dion J ( 27r ) 

pr ~ . . f d d p _ 



~ pf ~ r d p 

Qrep = -~ — , D rep (k) = 1- pr / j7T^9o(k - p)g (p) 
J-> rep J \^) 

(A.2) 

For the remaining index combinations we find that the pair symmetry is recovered and 
that Q(aa)(ab)(k) = Q( a b)(aa)(k) = Q(aa)(bc)(k) = Q \ bc ) (aa) (k) . The non trivial momentum 
dependence of these correlators can be understood since the q aa operator just checks for 
the presence of a spin, whereas the q a b is the squared magnetisation and obviously depends 
on the presence of a nearby spin. 
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The fu 

Qa = Q(aa 
QC = Q(aa 

Qp = Q(ab 
QQ = Q{ab 
Qr = Q(ab 



1 solution is, 

(aa) = Qb = Q(aa)(bb) = P 

(ab) = Qd = Q(aa)(bc) = PI + ^cQlc 



(ab) 
(ac) 



P + 2(S; on - ^rep)(Qlon ~ 3Q rep ) + 2Tj C dQc + ^RQlor}j 

~ (p + {Tjlon - ttrep + \){Qlon ~ SQrep) + ^CdQc + ^rQioti) 

)i„„ V ' 



D 



1 



(ac) = Qrep + ( P + (~ 2£ r ep + 1) (Qlon ~ 3Q rep ) + ^CdQc + ^RQlon ) 

(A.3) 

Where besides the longitudinal and replicon components some of the other E combinations 
are needed, 

d d p r Ad * 



/dp f dp 

j^y9o(k - p)gi(p) + \pq J j^y9o(k - p)(g - 2gi 



)(p) 



2(1 -D rep ) 

2"(1 — Di on ) 

X CD (k) = p(-2 + 2q + 3r) 
t R (k) = p(2 + 2q + 3r) 



d d p 
(2tt) c 



go{k-p)hip) 



d d p _ _ 

9o(k ~ p)go(p) + 3p(l - q - r) 



(27T) 



d d p 
(2tt) c 



+ ip(-l + 2 9 + r) J AlLg^k-pfafr) 



9o(k-p)gi(p) 



(AA) 



These final forms for each individual component are not especially illuminating, but it is 
worth noting that they have denominators only of the form -D; on , D rep , Df on and D rep Di on . 
The replicon denominator vanishes first, 



D rep (k = 0) = 1 - pr 



d d P ~2 ( \ 1 
9o(P) = 1- 



(27T) 



rgi 
9 



(A.5) 



So we recover the AT condition as the divergence of this correlator. 

The case in which q vanishes is particularly simple, D rep = Di on = D and only the P 
combination (and Q( aa )(bb) = p) remains non-vanishing, 

P 



Q(ab)(ab) (k) 



D 



(A.6) 



Similar techniques could be used to analyse the correlator in the RS but ferromagnetic 
phase and also for solutions with one step of RSB. 



30 



Appendix B. Stability of RS solutions 

We present the usual stability analysis based on the eigenvalues for small fluctuations 



about a RS solution along the lines of the work of de Alemeda and Thouless [19]. We 
therefore need the Hessian matrix H^ ab ^ cd y. 

<) ' l ' r '"' j\ oc (G- c \k)G^(k) + G-^k)G- c \k)) 5(k - k f ) - j^Vabcd (B.l) 



5G ab (k)6G cd (k') K acw d6W " dW 1 ' (27r) d 



For RS solutions the components (in the notation of [|19| where all indices are different) of 
the Hessian matrix are given by: 

A = H {aa)(aa) = 2 ( ^ - ^] 8(k-k')- 



C = H ( 





1 


2 






\9o 


2 ( 






\9o. 


H 


(ab)( 



(2Tl) d 



B = H (aa)(bb) = 2(^j 5(k - k') - 



(27T) 

2^ (^-^]s(k-k')--^- d n 

Wo J \9o 9 2 oJ (27v) d 



D = H (aaKbc) = H (bc){aa) = 2 (|) 5(k - k') - j£y^ ab (B.2) 
P ^^^((|) 2+ (^-|)l^- fc/) -(4° 



Q = H {ab)(ac) = - (!) (1 - 2 f) S(k - k>) - ^ bc 

R = -ff(ab)(cd) = 2 ( |^ j 5{k - k') - , P . d VL abcd 

Diagonalisation in replica space gives the replicon operator as the combination P — 2Q + R. 
It is convenient to introduce the positive weight function into the /c-space eigen- 
equation to obtain, 

m m ~ I idk ' m = K "m m (R3) 

Where r = (1 — 2Vt ab + Vt a bcd) is defined in ( |4.10| ). In the case where J f = 0, the eigenvalue 
is clearly positive, otherwise we integrate over k to obtain the condition, 

f d d k n r . 

pr j wr s °^v 31<1 (B ' 4) 
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The longitudinal mode is slightly more complicated because the diagonal terms 
A, B, C, D contribute. The equations for the eigenvalues are, 



(A - B)f + (A- B) - 2{C -D) h = X lon (fo + fi) 



2(C - D)f + 2 (C - D) + {P - 4Q + 3R) 



fl — Azon/i 



(B.5) 



With some rearrangement, and the same weight function, we find, 

(<7o - 2<7i 



9o 



— fo + fl — honifo + fl) 



f + p(-2 + 2q + 3r)g 2 J h(k') = X lon f 



;b.6) 



An inspection of these equations in the same way as for the replicon leads to the following 
condition for stability, 



p(-2 + 2q + 3r) 



d d k 
(2tt)° 



^0(^0 - 2<7i) < 1 



(B.7) 



For the diagonal RS solution in which q vanishes, r = 1, we find that both conditions 
( |B.4| ) and ( [B.7| ) are the same and we recover the same line fl4.7p for the spin glass transition. 
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